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Abstract 

A new formulation for the equation of motion of interacting dislocations is derived. 
From this solution it is shown that additional coupling forces, of kinetic and inertial 
origin, should be considered in Dislocation Dynamics (DD) simulations at high strain 
rates. A heuristic modification of this general equation of motion enables one to in- 
troduce retardation into inertial and elastic forces, in accordance with a progressive 
rearrangement of fields through wave propagation. The influence of the correspond- 
ing coupling terms and retardation effects are then illustrated in the case of dislocation 
dipolar interaction and coplanar annihilation. Finally, comparison is made between the 
modified equation of motion and a precise numerical solution based on the Peierls- 
Nabarro Galerkin method. Good agreement is found between the Peierls-Nabarro 
Galerkin method and the EoM including retardation effects for a dipolar interaction. 
For coplanar annihilation, it is demonstrated that an unexpected mechanism, involv- 
ing a complex interplay between the core of the dislocations and kinetics energies, 
allows a renucleation from the completely annihilated dislocations. A description of 
this phenomenon that could break the most favourable reaction between dislocations is 
proposed. 

keywords: dislocation interactions, elastic waves, plasticity of metals, inertial effect 

1 Introduction 

During the last decade, the understanding of crystal plasticity has been considerably 
improved with the rapid development of Dislocation Dynamics (DD) simulations. These 
simulations have indeed the capacity to quantitatively predict the behaviour of micro- 
metric samples from the modelling of the motion of discrete dislocations and interac- 
tions at the elementary scale 0] |2] [3] |4] [5) . Up to now, these simulations have been 
mainly devoted to low-strain rate deformations. In this case, dislocations motion is 
well described as a steady-state motion, the inertial forces being negligible. Recent 
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extensions toward more dynamic loadings (e.g. shock loadings) pointed out that iner- 
tial effects can be important, notably to overcome obstacles like other dislocations or 
defects @ [7] [8] [9] [10) . This is why, to extend the capacity of DD simulations to high- 
strain rates, works have been dedicated to the complex problem of dynamic equation 
of motion ifTTI [121 [131. 

Inertial effects for a single dislocation stem from the modification of the amount 
of energy, both elastic and kinetic, that follows change in the dislocation velocity. To 
balance such variations, supplementary work has to be done by the so-called inertial 
force. A simple estimation of this force relies on the hypothesis of steady-state stress 
and velocity fields around a dislocation |14|. However, rearrangement of the fields 
through wave emission has been shown to be critical in order to quantitatively describe 
the effect of inertia lTTTl[T3l . 

Wave propagation naturally leads to a retardation of the interaction between dis- 
locations that may have a very strong influence on shock loadings. In conventional 
DD simulations, the change of elastic forces due to dislocation motion is considered 
instantaneous, without being limited by the speed of the sound waves. Therefore, the 
dislocations moving behind a shock front artificially alter the stress in front of it. To 
avoid such unphysical propagation, retardation effects have to be considered in the 
dislocation-dislocation interaction forces. 

In this paper, we propose a modelling for inertia and retardation effects in the 
framework of the equations of motion developed for the Dislocation Dynamic method. 
Two Equations of Motion (EoM) for interacting dislocations are proposed. In a first 
section, the solution of instantaneously updated fields allows the definition of all the 
terms appearing in the EoM. In this first equation, the overall kinetic energy does not 
reduce to the sum of the kinetic energies of each isolated dislocation. This induces 
additional coupling terms between dislocations of kinetic and inertial nature, the im- 
portance of which will be discussed. The second section is devoted to a heuristic mod- 
ification of the first EoM in which retardation effects are included in the inertial and 
elastic interaction forces. In the third section, the results from both EoMs are compared 
to the results of a full-dynamic and more fundamental method, the Peierls-Nabarro 
Galerkin method. This method corresponds exactly to the theoretical framework used 
for the proposed EoM, excepted for the treatment of acoustic waves, now exactly re- 
solved. This comparison is concluded by a discussion on the influence of inertia, on 
retarded effects and on the coupling terms accounted for in a dynamic EoM, for a 
dipolar interaction and a coplanar annihilation. In the latter configuration, an unusual 
mechanism of renucleation from the annihilated dislocations is depicted. A discus- 
sion of this phenomenon is proposed in the last section and emphasizes the complex 
interaction between kinetic and core energies of the dislocations. 

2 Instant Equation of Motion 

Many mechanisms involved during dislocations interactions can be investigated with 
the simplistic problem of two attractive parallel dislocations of opposite signs (noted 
a and j3). In this model, fields around dislocations are supposed to be modified ev- 
erywhere in a time interval very short compared to the time needed by the acoustic 
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waves to propagate in the solid. Therefore, at each time, the fields are close to sta- 
tionary solutions and are only function of dislocation location (e.g. x a ) and velocity 
(e.g. v a ) 03] • The total energy of the system E is obtained from the overall velocity 
fields u — u a + vfi and stress fields a = & a + & 13 , obtained by summing up the 
contributions of the two dislocations : 

E = e aa + + 2e al3 + k aa + k w + 2k a0 (1) 

where e denotes an elastic energy and k stands for a kinetic one: 

e"* 3 = \ f a a {r-x a ,v a ):C- 1 -.^{r^x^v^dr 
2 Jn 

k * 3 = -p t u a {r-x a ,v a )-u p (r~x p ,v p )dr 
2 Jn 



(2) 
(3) 



where C is the stiffness tensor. Terms noted by a double superscript D m or D' 3 ' 3 are 
related to isolated dislocations whereas mixed ones represent the cost of the interac- 
tion. With a hypothesis of instantaneous updated fields, the energy is only function 
of locations x a and x@ and velocities v a and . In the case of two dislocations of 
opposite signs and of symmetrical trajectories, energy conservation leads to the EoM 
(given in the following for the dislocation a): 

dx a v a dv a 1 1 

Terms in the right-hand-side of equation (|4]i are proportional to acceleration and are of 
inertial nature whereas terms in the left-hand-side represent interaction forces. A brief 
description of each terms will be given now. We note F ; E = -2de a P/dx a the classical 
elastic interaction force and = —2dk a0 /dx a a kinetic interaction force (the sub- 
script T stands for "instantaneous"). Inertia is made of two terms, the first one (termed 
self-inertial, SI) F ; SI characterizes the inertia of a single and isolated dislocation, as 
already defined by Hirth et al. [14|: 

SI = ^(e- + fc ««) = 

where m(v) is termed the instant mass of a single dislocation. This mass is a complex 
function of velocity and becomes unbounded for v approaching the shear wave velocity 
due to the divergence of strain and velocity fields (given in the following for an edge 
dislocation): 



m(v) = m s , (^) [-8tl - 20^ + A^ 3 + 7 7 s + 257^ - H7 S " 3 + 3 7s 

(6) 

with the mass of a screw dislocation at rest m Si o = ^ > depending on the 

1 /2 

shear modulus p and on the shear wave speed cs and 7l.s = (1 _1,2 / c ls) > w i tn me 
longitudinal wave speed cl. The instant mass depends on the size R of the domain in 
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which the strain and velocity fields are supposed to be adapted to the present dislocation 
velocity. The parameter ro is a cut-off radius usually chosen to be equal to b. 

The second inertial term F; 11 = 2(v a /v a ) d(e al3 + k al3 )/dv a , represents an"inter- 
inertial" (II) force. It can be noted that the equation of motion (termed in the 
following as the "instant" EoM), which now reads 

Fj E + F; K = if 1 + F/ 1 (7) 

contains two coupling terms (i^ K and F ; n ), usually not considered in studies of inertial 
effects for high velocity dislocations |[T6l[T7l[T0l . However, the influence of these terms 
has not been shown to be negligible. In particular, when two dislocations superimpose 
(for example when a junction is created), the overall energy is not reduced to the sum 
of the individual energies. Coupling energies e"' 9 and k a/3 , from which the forces Fj K 
and i-f 1 are derived, can not been neglected in general and could a priori play a role. 

An estimation of the kinetic and elastic interaction energies is obtained numerically 
for two straight and parallel dislocations of opposite signs and velocities. By noting that 
with the considered symmetries x a — —xp, equations (|2]i and (0) reduce to a simple 
convolution product performed with a fast Fourier transform. Relativistic stationary 
fields for stress, strain, velocity of a finite core dislocation are considered 1 15 1. 

The hypothesis of instantaneously updated fields is however a quite strong assump- 
tion since any changes in the velocity of the dislocation cannot be propagated more 
rapidly than the shear or longitudinal wave celerity. This is particularly true for high 
strain rate loadings in which a significant dislocation motion can occur during this 
propagation time. Hence, in the following section we construct a modified EoM in 
which retarded effects are now considered. 



3 Retarded Equation of Motion 

In this section, we propose a heuristic modification of the elastic and self-inertial forces 
in which retarded effects are included. The retarded forces F r n and F r K are far much 
difficult to derive than F; 11 and F ; K because of the hypothesis of unsteady velocity and 
stress fields. For this reason no numerical estimation of these terms is given here. 

3.1 Retarded inertial force 

To get rid of the assumption of the stationarity made in the first model, we use a self- 
inertial force that takes into account emission and propagation of waves accompanying 
changes of velocity of the dislocation fl3l . This solution is constructed by using the 
retarded self-inertial force (F r SI ) produced at time t by an instantaneous velocity jump 
from to v at time r < t fTZl : 

F^(t-r,v) = ^- . (8) 

t — T 

with g a function that depends on dislocation character. The work done by this trailing 
force balances the increase of total energy due to the progressive updating of strain 
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and of velocity fields from the solution at v = to the one at v > 0. The retarded 
inertial force for any function v(t) is constructed by summing all the contributions of 
the trailing forces SF — [dF(t — T,v(T))/dv(T)]Sv(r) due to elementary velocity 
jumps Sv at t = r which we assume to be a reasonable approximation of the trailing 
force due to a jump at r from v to v + Sv: 



This expression of the self-inertial force is singular at r = t due to the assumption of 
point dislocation done in the original work 111211 . A regularization of the time-kernel 
has been proposed by Pillon et al. [13] to account for a core-size for the dislocation 
and the retarded-inertial force now reads: 



where to = Co/ c s and Qq the dislocation core width at rest. In [13] it has been un- 
derlined that this expression leads to extremely small inertial force (compared to the 
one given in equation ©) when one focuses on short times scales, and has a complex 
dependence with respect to time, according to the fact that the zone experiencing an 
energy evolution is progressively expanding from the vicinity of the dislocation to the 
whole solid. This non-local in time force represents the interaction of a dislocation 
with its own past motion. For interacting dislocations, the construction of a retarded 
elastic interaction force is now proposed. 

3.2 Retarded elastic force 

The retarded elastic interaction force P t E is due to stress around moving dislocations 
that progressively rearrange by wave propagation. We propose a construction of re- 
tarded elastic forces between moving dislocations. In the case of a single dislocation 
that jumps from v to v + Sv at instant t = r, the stress field at a distance \p — x(t)\ 
from the dislocation has been modified by the velocity jump Sv, provided that acoustic 
waves have been propagated up to this point, which is the case of any point p satisfying 
the condition \p — x(t)\ < c(t — r). We suppose that the stress at this point is equal to 
the stationary solution at velocity v + Sv (this is the case of the point "Pi" in the fig- 
ure[lj-left). Conversely, the field at distance \p— x(t) \ > c(t—r) from the source of the 
waves is not modified by the jump. It is therefore the field of a dislocation in stationary 
motion that ignores the velocity jump, that is to say with a velocity v and located at 
x(t) + (t — t) ■ v, different from the present position x(t) = x(t) + (t — r) • (v + Sv). 
This second situation is the one of the point "P2" of the figure [T]-left. Therefore, two 
stationary solutions have to be considered, depending on the relative location of the 
examination point (in "Pi" or in "P2", see figure [T]-left) where the stress is evaluated 
and of the dislocation trajectory. 

One can notice that the stress at P2 is produced by a virtual dislocation that, from 
the time r has conserved the same stationary velocity v(j) (i.e. Sv = 0) up to the 
present time t. In the case of a constant velocity the virtual dislocation is superposed to 




(9) 




(10) 
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Virtual stationary dislocation 




Space x ( % ) Space 



Figure 1 : Definition of the retarded elastic force. The solid curve represents the tra- 
jectory in a space-time diagram and the dashed-dots lines, the trajectory of the virtual 
dislocation which would have kept a stationary velocity. On the left-hand side is shown 
the delay induced by a jump from a stationary velocity v to another velocity v + 5v and 
on the right-hand side, a general trajectory is illustrated. 



the real one and corresponds to the standard way to calculate interaction forces between 
moving dislocations lfl0l[T6l . 

For any dislocation motion, we suppose that a stationary solution is achieved at 
each time step but is visible only after the propagation of the waves and before any 
new change in the velocity, which will lead to a new stationary solution. Thus, the 
stationary solution to be considered at time t and at the measurement point p is given by 
the element of trajectory in the past time r for which the relation \p — x(t) = c(t — r) 
is verified, as depicted in figure QJright by the gray zone. 

In the case of two interacting dislocations of circulation £, the Peach-Kohler inter- 
action forces (F^> °^ a (t) coming from /3 to a and conversely F®' a ~ Jf P(t) from a to 
/3)is 

Ff-^it) = {a [Ap al3 , v {t )] -b a }A£ a (11) 

with A is the cross product, Ap a/3 = x a (t) — [x l3 (T l3 )+v^(T l3 )-(t — T 13 )] represents the 
relative location between the dislocation a and the virtual dislocation j3 and <r [Ap, v] 
the stress tensor at Ap of a dislocation with a velocity v. The past time r " at which the 
stationary solution is considered is given by the intersection of the trajectory x^{t^) 
with the space-time cone \x a (t) — x^{t^)\ = c(t — r* 3 ). 
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The force F®' a ^ r ^{t) is obtained by a similar way but may result in time r Q and 
distance Ap a/3 that is distinct from r* 3 and Ap* 3 ". This is different to the standard 
definition of interacting forces for which no propagation time are considered (that 
is = r a = t), leading to an equality of the distances between dislocations (i.e. 
Ap a P = Ap^). 

In this simplified construction we use stationary solutions, which amount to the 
neglect of transient waves accompanying any changes in velocity. Removing this hy- 
pothesis is made possible by considering exact solutions, as for example in the method 
based on Green functions proposed by Mura [18| for expanding loops. These exact 
solutions involve however an additional temporal integration of a prohibitive computa- 
tional cost for DD simulations. 

4 Applications 

The EoMs of the previous section have been implemented in a two dimensional (point 
like) DD code for validation. Such validation is made by comparing the results of 
the DD simulation and those of a more fundamental nature called Peierls-Nabarro 
Galerkin (PNG) simulation |fl9l |20l . Indeed, the latter technique allows for a full- 
dynamic description of dislocations interactions and has the advantage to share the 
same set of hypothesis than the ones used for the definition of the EoM, namely, an 
isotropic elasticity, a continuous representation of the solid and a simplified micro- 
scopic viscosity. 

A first simulation (a dipolar interaction), dedicated to the analysis of retardation 
effects is proposed in which the coupling terms and F^ 1 are expected to play a 
minor role. A second simulation (a coplanar annihilation), which can be seen as a two- 
dimensional substitute for a more general study on junctions magnifies these coupling 
terms, with however a non-negligible influence of retardation mechanisms. 

For each configuration, the two dislocations are supposed to be parallel edge dis- 
locations of opposite signs. The line direction of the two dislocations are oriented 
along the z-direction and they have opposite Burgers vectors of magnitude b. They 
are separated along the y-direction by a distance h and are restricted to glide in the 
x-direction. A viscosity is introduced on the glide plane for the PNG simulations and 
permits the existence a finite dislocations velocity for constant applied stress. This 
"microscopic" viscosity brings to a "mesoscopic" dislocation viscosity, following the 
model presented by Rosakis lETI and is implemented in the DD technique. Since the 
aim of this paper is to analyse and to model retardation and inertial effects, whatever 
the model of viscosity considered, no attempt was made to introduce complex dissipa- 
tion phenomena emerging from the interaction of the dislocation core with the atomic 
lattice [22 23 24] . In the following subsection, we briefly describe the PNG methods 
and some new improvements developed to provide comparisons with DD simulations. 

4.1 Reference PNG simulations 

The Peierls-Nabarro Galerkin (PNG) method [19 20 1 is a generalization of the Peierls- 
Nabarro concept in which the displacement fields are represented by an element-free 
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Galerkin method, close conceptually to the finite elements method. The non-linear be- 
haviour is introduced by allowing a displacement jump 77 along the glide plane at the 
cost of an additional energy 7 lsf (the inelastic stacking fault energy) that is deduced 
from the 7— surface [19|. Incorporation of kinetic energy allows for acoustic waves, 
which are essential for instationary dislocations motion (see, for example the modelling 
of accelerated dislocations iTHl ). This method has shown to reproduce very well the 
analytical solutions for a stationary-moving dislocation, even in the high velocity (rel- 
ativistic) regime and has been used by Pillon et al. |[T3ll to check the EoM defined by 
equation ( fTOt for a single dislocation. 

All the DD simulations are done in an unbounded domain f2, naturally avoiding 
dislocations images [25 j. To prevent from dislocation images in the corresponding 
PNG simulation, the displacement fields on the boundary dVl is given by the convolu- 
tion product of displacement field of a point dislocation in stationary motion with the 
dislocation density Vry(s) 

u imp (r dn ) = Vri(s)u stat (r dn - 8, v(s)) ds , (12) 

where L is the glide plane. The displacement field u stat (r, v) is the exact relativis- 
tic displacement field [15] generated by stationary dislocation moving at the velocity 
v(s). The straightforward choice for v in equation dT2b should be the velocity of each 
infinitesimal dislocation v = fj/Vrj. Nevertheless, this definition induces a noisy mea- 
sure of v(s) which is transmitted to the boundary conditions and may generate acous- 
tic waves. The interaction of these waves with the dislocation modifies the velocity 
v(s) and can bring on increasing oscillations. To get rid of this possibly resonating be- 
haviour, the velocity v(s) is replaced by an average along the glide plane of the velocity 
of each infinitesimal dislocations v = (7//V77). This averaged velocity v is then filtered 
in time with a first-order filter v + TfV Rlt — v and is used to define the imposed dis- 
placement u lmp . The characteristic time T{ is set to the time needed to accelerate the 
dislocation up to stationary motion. 

To be comparable with the DD simulations, PNG simulations of the interaction 
must be done with dislocations that are close to stationarity. To clean all the instationary 
information due to the initial acceleration of the dislocations, a body force proportional 
to the difference between a stationary velocity field tt lmp and the present one u is 
applied 

div (<t) — pu= — — u inip (r n )] , (13) 

where t c is a characteristic time and tt lmp (rn) the field defined in equation ( fT2b in the 
whole volume f2. The convergence time towards the stationary solution is of the order 
of r c . 

This convergence mechanism is set to zero (ie r c — > 00) at the beginning of the in- 
teraction so that the instationary solution is no longer altered. Thanks to this procedure, 
displacement and velocity fields are in accordance to stationary motion. The compari- 
son with DD simulations presented in the following begins after this initialization step. 
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4.2 Dipolar interaction 

The coupling terms F ; K and represent the difference of kinetic and inertial energy 
between interacting dislocations and the same dislocations considered as isolated ones. 
Preliminary simulations have shown these terms to decrease rapidly when the distance 
between dislocation increases. A dipolar interaction with a distance between gliding 
planes of several Burgers vectors is therefore a good configuration to test retardation 
effects, with a weak influence of coupling terms. Consequently, the kinetic interaction 
force and the inter-inertial force can be neglected in the dipolar interaction problem. 

The initial trajectories and boundary conditions used in the DD and PNG simula- 
tions are the same with a steady state motion at t = 0. The two dislocations move 
symmetrically and at each time, the retarded elastic force is calculated by searching in 
the past the element of trajectory defining the location and velocity of the virtual dis- 
location, as already described in section[3] The distance h in the y-direction has been 
fixed in PNG and DD simulations to 86 which corresponds to an average value of min- 
imum heights for dipole observed experimentally ||26l . The cut off radius R = 500/irn 
(or 25006) in equation © corresponds to a typical dislocation density of 10 12 /m 2 . 
The influence of retarded elastic force, and of retarded inertial forces are analyzed by 
turning on and off Fp/F^ and F? l /F^ 1 . In F®, the velocity c is supposed to be the 
velocity of shear waves eg, in accordance with observations of Pillon et al. [ 13 1 where 
eg brings the main contribution to self-inertia. 

Figure |2] shows results obtained by the PNG technique and the DD technique, the 
retarded aspect for the F E and F sl being turned on and off. The simulations are done 
for a constant applied stress cr a . During the short-range interaction, DD simulations 
show no influence of the retarded interaction force. This is consistent with a propa- 
gation time of about 86/c (or in dimensioned time « 0.3 ps) short compared to the 
characteristic time of the interaction (w 506/c, or in dimensioned time ss 1.7 ps), in- 
ducing negligible retardation effect. This similarity progressively vanishes when the 
interaction distance increases, and the instant elastic force eventually gives a stable 
configuration (v = 0) whereas the retarded force predicts a complete separation of the 
two dislocations. The corresponding PNG simulations cannot be achieved up to this 
time but clearly follows the retarded interaction simulation. 

The difference between EoM with retarded self-inertial term and EoM with an in- 
stant mass is more contrasted during the short-range interaction, the better match with 
PNG simulations being obtained by the EoM with retarded self-inertial term. This dis- 
crepancy is due to an overestimation of the characteristic size of the zone playing a role 
in the instant mass. In the instant EoM, this size has been set to R = 25006. However, 
most of the variations of the velocity take place in a time range of less than 50 6/c s , 
which limits the zone contributing to inertial effect to 50 x 6 (see figure [J), leading to 
a strong overestimation of the inertial effect by the instant self-inertial force. Finally, 
the critical stress that breaks the dipole is found to be 0.66%/i without inertial effect 
and reduced to 0.42%/i by using the retarded EoM whereas the instant mass gives only 
0.37%^. 

Therefore, retarded inertia play a part mostly when two dislocations cross each 
other (or, in a more general case when a junction forms) and cannot be modelled accu- 
rately with an instant mass. The dipole formation occurs thus easier than expected with 
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Figure 2: (Colour online). Velocities vs time for a dislocation interacting with another 
one in a dipole of edge dislocations. We first consider EoM without retardation (F ; E = 
-FJ ), then we add the retarded inertial force (F ; E = F^ 1 ) and finally we use fully 
retarded EoM (F t E = if 1 ). 
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an instant mass EoM, which promotes the crossing over between dislocations by an 
overestimation of their kinetic energies. Conversely, the effect of retarded elastic force 
influences mainly the long range interactions, and intensifies as the distance between 
dislocations increases. This can potentially have a strong influence for high strain rate 
since the adjustment of the dislocations location is only perceived in their local envi- 
ronment. This is especially relevant for shock loadings in which the stress modification 
induced by dislocations motion should be confined behind the shock front. 

4.3 Annihilation 

The use of a coplanar annihilation in place of the dipolar interaction results in an im- 
portant increase of the relative influence of the coupling forces F K and F u . A stiffer 
variation of the elastic forces (notably when the dislocations are superposed), leads, fur- 
thermore, to a significant sensitivity to retardation effects. Surprisingly, this reaction, 
which can be considered as the strongest possible one Q, is shown in the following to 
be also breakable by inertial effects. This configuration is therefore a more severe test 
of the proposed EoM than the dipolar one, with an additional difficulty coming from 
the impossibility to separate the influence of each of the forces, coupled and retarded. 

An additional interaction force, specific to coplanar annihilation, comes from the 
possibility to superpose the dislocation, by the way modifying the overall energy stored 
into the dislocation core. Indeed, in the framework of the Peierls-Nabarro Galerkin 
method, the energy of a single isolated dislocation contains a core energy J 7 lsf [r)(x)] dx 
with 7 lsf a surface potential and rj(x) the displacement jump along the glide plane. This 
energy is constant for dislocation with a fixed core width and its influence is usually 
neglected. However, this energy completely vanishes when two dislocations of oppo- 
site sign superimpose. The corresponding potential energy is obtained by summing the 
displacement jump of the two dislocations : 



where 7 is the interplanar potential used in the PNG method and where x a and x" are 
the dislocation locations. A rough estimation of this interaction force consists in taking 
for rj(x) the quasi-static solution of the dislocation displacement jump. An additional 
force F 1 = —d r y a ^/dx ci derives from this potential, and appears in the left-hand side 
of equation (0. Therefore, both influence of coupling forces, retardation effects and 
core energy could have an influence in this configuration and will be measured. 

We extend the method used for dipoles to coplanar annihilation by fixing the dis- 
tance between glide planes to h = and by introducing a finite core size to avoid 
singularity in the elastic force when the two dislocations meet. Contrary to the dipolar 
interaction, we measure the minimal applied stress above which dislocations renucleate 
after annihilation. The overall result depends potentially on coupling terms, retarda- 
tion effects and core energy. The influence of the coupling forces F ; K and F^ 1 is tested 
by using the "instant" EoM and by turning them on and off. The role of retardation 
effects is estimated by comparing the retarded EoM (in which no coupling forces are 
known) and the instant one, in which F ; K and i*] 11 are switched off. The influence of 
core energy F 1 will be studied in the two equations of motion. 




(14) 
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Results obtained in each cases are summarized in figure [3] The F 1 is found to 
modify only the retarded EoM. This is consistent with an evolution of this force only 
when the dislocation cores are in contact. This force is therefore a very brief signal 
that is felt only with the retarded EoM which is known to predict low inertia for high 
frequency loadings [ 13"]. 
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Figure 3: Summary of results obtained for the coplanar annihilation. 



Contrary to the force F 1 ', the introduction of F and F in the "instant" EoM de- 
creases slightly the critical applied stress. Indeed, the energy k a/3 is positive and has its 
maximum when e al3 is minimum. The kinetic interaction force partially compensates 
the elastic interaction forces, which in turn decreases the critical stress. The force F u 
has approximately the same influence as F K . Actually, for low velocity dislocation 
the derivative of e aa + e a/3 with respect to the velocity is very low (the elastic field 
weakly depends on velocity) and only the energy k aa + k a/3 intervenes in the mass. 
For the same reason as before, k aa + k a ^ is intensifying during annihilation, by the 
way increasing the inertial energy which helps the crossing of the dislocations. 

The critical stresses given by the proposed EoMs are however far from the results 
of PNG simulations that predict a crossing for applied stresses of 2.7%/i but a complete 
annihilation for 2.4%/j. An explanation of this discrepancy is now proposed. 



5 Discussion 

Quasi-static reactions between dislocations are known to be correctly described as lines 
in elastic interaction without any core contribution 11271 12]. It is striking to notice that, 
in opposition to the quasi-static case, dynamic annihilation cannot be quantitatively 
modelled by the proposed EoM, even by taking care of wave propagation through re- 
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tarded mechanisms. The dislocations energy is therefore damped down by another 
mechanism that cannot be represented in terms of superimposition of two dislocations. 

The difference between analytical modelling and PNG simulations can be explained 
by analyzing the role played by the interplanar potential 7 Isf during the reaction (fig- 
ures |4] and |5). In figure |4]-a, the two dislocations are getting closer at a velocity of 
» 0.83c s and eventually superpose (figure[4]-b). At this time, no dislocation is present 
(i.e. 7] w everywhere) and the energy is mainly of kinetic nature. The available 
amount of energy results in the creation of two dislocations of magnitude half the ini- 
tial Burgers vector in figure [4]-c. This two "partial" dislocations, separated by more 
than 4 6, are thus accumulating, as a staking fault, an important amount of potential 
energy. This state is very different from the expected "ideal" entire dislocations (half 
of each dislocations is missing) and changes dramatically the displacement fields that 
induces an intense stress wave (figure 0]-c). This state is however unstable since the 
7 lsf does not have any local minimum at 6/2 (see figure|5]): the field rj has to go through 
this potential barrier to achieve rj = 6 or to turn back tor] = 0. 

For an applied stress of 2.4%//, kinetic energy is not high enough to allow the 
dislocation to overcome the maximum of the inelastic stacking fault energy at rj = 6/2. 
This is illustrated on figure [5] where the inelastic staking fault is plotted with respect 
to the position on the slip plane and to the value of the displacement jump rj. Since 
the available kinetic energy is not sufficient to overcome the barrier, the displacement 
field rj minimizes the potential energy by decreasing to a homogeneous value of 77 = 
along the slip plane. To compensate this decrease of the inelastic stacking fault, a 
second wave is emitted in the solid (see figure @]-right-e). Conversely, for an applied 
stress of 2.7%/x, kinetic energy is large enough to overcome the potential maximum 
for rj = 6/2 and to create two dislocations (with a displacement jump ranging from 
77 = to 77 = —6). The sudden change from two half dislocations to two entire ones 
induces the emission of a wave in the solid (figure|4]-left-d and -e). The "renucleated" 
dislocations moves with a velocity of w 0.78c s . 

Therefore, most of the available energy is dissipated in a "two steps mechanism". 
A first step consists in storing an important part of the kinetic energy into staking 
fault energy while another part is lost by an intense acoustic emission. The second 
step results in a complete relief of the staking fault energy, leading one more time 
to an important acoustic wave emission. This scheme, very different from a simple 
superposition of dislocations, explains the discrepancy between PNG and the EoMs. 
Such phenomenon could also a priori intervene during other contact reactions like 
junctions formation, as far as junction breaking due to inertia can be invoked. 

6 Concluding remarks 

An equation of motion for interacting dislocations is proposed by using two descrip- 
tions for stress and velocity fields, denoted "instant" and "retarded". A comprehensive 
study of the forces acting during dislocations interaction shows that in addition to the 
usual elastic and inertial terms, a kinetic interaction force and an inter-inertial force 
should be considered in EoMs for fast interacting dislocations. 

We show that inter-inertial force and inter-kinetic force does not play a significant 



13 



t=3.30ps 



t=3.37ps 



33%H 

* * I' 



-20 -10 10 20 



Space/b 



- - 



-20 -10 10 20 



t=3.65ps 



-20 -10 10 20 



Space/b 



t=3.74ps 



Space/b 



I 



Q 



-20 -10 10 20 

t=4.04ps 



Space/b 



-20 -10 10 20 



Space/b 



Space/ b 



-20 -10 10 20 



Space/ b 



-20 -10 10 20 -20 -10 10 20 

=3.91 ps , 



Space/ b 



t=3.93ps 




n 


9 


d 


. Space/b 


-20 


-10 





10 


20 


t=4.10ps 


( 




1 
1 


e 


. Space/b 



-20 -10 10 20 



Figure 4: (Color online). (Left) Snapshots of PNG simulations of a coplanar interaction 
with 2.7%/i as applied stress. Pictures represent the equivalent stress in the whole space 
with in addition the corresponding displacement jump on the glide plane. (Right) The 
same simulation with 2.4%// as applied stress. 



role during short distance interaction like formation of dipoles, as far as an "instant" 
EoM can be considered. The retarded effects introduced in the inertial terms become 
important for contact reactions like annihilation and possibly formation of junctions. 
For instance, the "retarded" EoM leads to critical stresses for dislocation annihilation 
two times higher than the "instant" EoM. 

In addition, we show that the elastic retarded force modifies long range interactions 
and therefore is essential to model shock loadings. In such conditions, moving disloca- 
tions will concentrate stress only behind the shock front and as a consequence nucleate 
original plastic features. 

The comparisons between DD and PNG simulations shows however that "instant" 
as well as "retarded" EoMs, are both failing to reproduce quantitatively the inertial 
effects observed at the limit case of coplanar annihilation. In such case, an original 
mechanism of energy accumulation into the interplanar potential is proven to be the 
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Figure 5: (Color online). Evolution of the inelastic stacking potential 7 lsf (77) along the 
displacement jump r\ and along the slip plane for different times. The two configu- 
rations corresponding to an applied stress of a = 2.4% fj, and a = 2.7% fj, have been 
represented. 



reason for DD simulations deficiency. More generally, we show that inertial effects can 
strongly influence contact reactions. As an example, two dislocations with opposite 
Burgers vector can completely annihilate and renucleate as a result of inertia. From 
this observation, related to the most energetically favourable dislocation-dislocation 
reaction, one can conclude that inertia may be determinant in many strain-hardening 
mechanisms involved during high-strain rate loadings. 
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